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EXPERIMENTAL  OBSERVATION  OF  THE  STRUCTURE  OF 
SHOCK  WAVES  IN  DUSTY  GAS^ 

Yu  Hongru,  Lin  Jianmin,  Yuan  Shengxue ,  and  Li  Zhongfa 
Institute  of  Mechanics,  Chinese  Academy  of  Sciences 

Abstract 

The  structure  of  shock  waves  in  dusty  gas  was  studied  by  shock 
tube  experiments.  Pressure,  shock  wave  velocity  and  dust  concentra¬ 
tion  were  measured  under  the  condition  that  the  dusty  loading  ratio 
was  less  then  0.4  and  the  shock  Mach  number  was  less  than  1.5.  The 
frozen  and  equilibrium  pressures  behind  the  shock  wave  can  be  well 
described  by  the  RankinerHugoniot  relation.  But  the  relaxation  length 
was  significantly  shorter  than  that  calculated  from  the  drag  of  a 
sphere  at  steady  motion.  And  the  disperse  shock  waves  in  dusty  gas 
were  observed.  It  was  found  that  the  presence  of  dust  decreases  the 
pressure  fluctuation  in  the  flow.  The  propagation  velocity  of  the 
shock  wave  agrees  with  the  result  of  the  model  analysis  on  an 
‘effective  gas*. 


This  paper  was  recommended  by  editor  Ling  Dong  Ji  and  was  received 
on  July  1,  1983. 

This  paper  is  written  for  the  15th  anniversary  of  the  death  of 
our  adviser  Professor  Guo  Yong  Hui. 


The  List  of  Symbols 


Definition 

Sound  velocity 
Dusty  gas  constant 
Particle  cross-section  area 
Time 

Heat  capacity  of  particle  material 

Gas  temperature 

Particle  temperature 

Gas  heat  capacity  at  constant  pres 

Gas  velocity 

Coefficient  of  resistance 
Particle  velocity 
Particle  diameter 
Voltage 

Coefficient  of  light  attenuation 
Distance 

Coefficient  of  heat  conductivity 
=x  le 

Intensity  of  transmitted  light 
Sound  velocity  of  dusty  gas 
Coefficient  in  formula  (2) 

-  *tT, 

Heat  conductivity 
Ratio  of  gas  heat  capacities 
Thickness  of  particle  cloud 
Specific  heat  ratio  of  dusty  gas 
Mass  flux  of  gas 


=C/C 

P 

Mass  flux  of  particle  cloud 
Gas  density 

Mach  number  of  shock  wave 
Density  of  dusty  material 
Number  density  of  dust 
Concentration  of  dust 


Nu  =hD/K  Nusset  number 

».  Velocity  relaxation  time 

r  Pressure 

It  Coefficient  of  powderiness  of  gas 

R  Gas  constant 

I?  =mp/m 

Subscript 

0  Wave  front  or  initial  value 

d  Dark  value 

f  Frozen  value 

c  Average  value 


1.  Introduction 

Intensive  theoretical  studies  and  numerical  analysis  of  the 
structure  of  shock  waves  in  dusty  gas  have  been  made  by  many  people 
[1-14].  Their  work  was  carried  out  with  some  idealized  assumptions. 
Among  them  the  most  basic  assumption  was:  the  particles  distribute 
uniformly  in  a  perfect  gas  medium  and  the  particle  cloud  was  treated 
as  continuous.  For  energy  and  momentum  exchanges  between  the  two 
phases,  only  the  drag  of  an  isolated  rigid  sphere  at  steady  motion 
and  heat  conduction  by  convection  are  considered.  Besides  that,  the 
results  of  experimental  observation  of  the  structure  of  shock  waves 
of  dusty  gas  are  rare  except  the  experiments  for  the  measurement  of 
the  effective  drag  coefficient  of  the  particle  cloud. 

In  this  paper  we  report  the  observed  results  of  such  structure 
of  shock  waves  in  a  shock  wave  tube,  and  the  driving  characteristics 
of  a  dusty  gas  shock  wave  tube. 

2.  Description  of  Experiment 

The  shock  wave  tube  for  the  experiment  was  installed  perpendic¬ 
ularly  with  driving  section  in  the  upper  position.  The  driving 
section  is  1.5  meters  long  with  an  inner  diameter  of  84  mm;  the  driven 
section  is  5  m  long  with  an  inner  diameter  of  185  mm.  The  reason  for 
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using  such  big  tube  is  to  decrease  the  disturbance  of  the  tube  wall 
to  the  minimum.  The  observation  area  is  located  4  m  below  the 
isolating  diaphragm. 

The  dust  material  is  aluminium  oxide,  with  a  density  of  3. 9-4.0 
g/cm^,  and  a  heat  capacity  of  0.73  KJ/kg-K.  The  shape  of  the 
particles  is  close  to  sphere  but  the  surface  is  not  smooth.  The 
dusty  particles  were  jetted  by  attraction  into  the  shock  wave  tube  in 
the  vicinity  of  the  isolating  diaphragm.  The  uniformity  of  the 
dusty  particle  distribution  and  the  initial  disturbance  have  important 
influence  on  the  experiment.  The  dusty  particles  were  jetted  upwards, 
then  the  particles  were  reflected  from  the  isolating  diaphragm  and 
fell  down.  Thus  the  movement  of  the  dusty  particles  was  basically  a 
free  fall  motion  and  the  distribution  was  quite  uniform.  Only  a 
small  amount  of  gas  is  needed  for  the  attractive  jet  and  the  initial 
disturbance  is  very  small.  All  the  gas  used  in  the  experiment  was 
nitrogen.  The  front  of  the  shock  wave  was  atmospheric  pressure. 

Besides  the  routine  static  measurement  of  the  initial  state, 
the  propagation  velocity  of  shock  waves  along  the  tube  axis,  the 
pressure,  and  the  density  of  dusty  particles  were  measured  in  the 
experimental  process.  The  light  attenuation  method  was  used  for  the 
measurement  of  dust  density.  When  the  sizes  of  the  particles  are 
uniform  or  the  size  distribution  is  narrow,  this  is  a  simple  and 
reliable  method.  We  chose  a  low  power  He-Ne  laser  as  light  source. 

The  number  density  of  dust  can  be  calculated  according  to  Bouger's 
law  by  measuring  the  intensity  of  transmitted  light  which  passed  the 
particle  cloud  [13]: 

^—exp—(n,A,lE).  (1) 

/t 


The  transmitted  light  was  measured  by  a  photo  diode.  The 

output  voltage  is  proportional  to  the  intensity  of  incident  light 
intensity.  Therefore  formula  (1)  can  be  modified  to  (2): 


<t(0  —  V-n 


v,-v^ 
no  -  Vi 


The  coefficient  k  depends  on  the  size  and  shape  of  dust 
particles,  the  thickness  of  the  particle  cloud  and  the  optical  loss 
in  the  optical  route.  We  determined  k  value  by  calibration  measure¬ 
ment  of  different  types  of  dusty  particles.  The  most  serious  problem 
in  the  measurement  of  the  density  of  the  dust  density  was  the 
adsorption  of  particles  on  the  observation  window,  which  decreases  the 
intensity  of  the  transmitted  light.  When  both  the  density  and  the 
velocity  are  low,  the  adsorption  will  be  low.  Such  measurement  of 
density  in  this  work  was  used  to  monitor  the  axial  uniformity  of 
dust,  and  to  offer  the  initial  values  of  the  dust  density. 

All  the  dynamic  measurements  were  recorded  by  instant  wave  form 
recorders,  and  the  data  or  graphs  were  output  after  the  signal 
analysis  instruments. 

3.  The  Positive  Shock  Wave  of  Dusty  Gas 


If  the  density  of  dusty  particles  is  not  high,  the  volume 
occupied  by  the  dusty  particles  is  negligible.  Then  the  equations 
of  the  conservations  of  mass,  momentum,  and  energy  are  as  flows: 


m  — 

m,  —  ijm  —  <7«,  — 
mu  +  i\mUf  +  P  —  (1  -p  n)mu,  -f-  P, 

C,T +  +  ,)  -p  C/J  +1,5)7. 


(3) 

(O 

(5) 

(6) 


and  plus  the  equation  of  state  of  ideal  gas: 


P-  pFT  . 


(7) 


In  order  to  solve  the  seven  variables  of  u,p,T,P,Uf,o  and  7, 
another  two  equations  are  needed. 


Generally  the  shock  wave  can  be  divided  into  three  regions  and 
disposed  separately: 


In 

density 


the  frozen  region  close  to  the  wave  front,  because  the  dust 
is  10^  times  higher  than  the  gas  density,  and  the  wave  front 


is  very  thin,  then  the  time  for  particles  to  pass  through  it  is  very 
short;  therefore  the  particle  velocity  and  temperature  can  be  consid¬ 
ered  frozen,  that  is: 


( 8 ) 

-  T.  ( 9 ) 

Substituting  these  two  equations  into  equation  (5)  and  (6) ,  all 

the  terms  related  to  the  dusty  particles  are  calceled  and  the  equations 
become  shock  wave  equations  of  pure  gas.  Then  the  gas  parameters  can 
be  obtained  by  solving  these  equations. 

In  the  equilibrium  region  below  the  wave  front,  the  equilibrium 
between  the  two  phases  is  attained  via  the  exchanges  of  momentum 
and  energy;  therefore  we  have: 

(10) 

T^-T,  ,  (11) 

Substitute  these  equations  into  (5)  and  (6)  and  reorganize 
these  equations.  Then  we  can  know  that  the  differences  from  the 
pure  gas  shock  wave  equations  are  as  follows:  Use  the  equilibrium 
sound  velocity  a,  —  [  +~j~]  replace  a^,  and  use  the  ratio  of 

heat  capacities  of  dusty  gas  r  — (l+  jj5)r/(l  +  rijs)  to  replace  r. 

Therefore  the  parameters  of  the  equilibrium  region  also  can  be 
obtained  from  the  pure  gas  shock  wave  equations. 


For  the  pressure  in  the  frozen  and  equilibrium  regions,  the 
following  relations  exist: 


El  -  _  y- 1 

P,  r+  I  \  aj  r  +  1 

ft  _  (HiV  _  r-  1 

Po  r  +  1  Vo,/  r  -t  1  * 


(12) 

(13) 


For  the  gas  velocities  in  the  frozen  and  equilibrium  regions, 
the  following  relations  exist; 


"ft 


"ft 


2rR 
r  +  1 


T,-l- 


T  —  1 

r  +  1 


2TR 

r  + 1 


Tt-h 


r  —  1 
r+  1 


"i 


(14) 

(15) 
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Because  “»%  T  and  R  are  always  smaller  than  ao-.  r  and  R,  therefore 
Pc  > /*/»  “r  <  “/.  According  to  the  equation  of  mass  conservation,  we  have 
Pr>  pt.  As  for  the  temperature  in  these  two  regions,  the  difference 
is  not  large,  and  which  is  higher  depends  on  the  particle  size  and 
gas  heat  capacity. 


For  the  two  lacking  equations  for  the  solution  of  the  relaxation 
region,  we  can  use  the  momentum  and  energy  exchange  equations  between 
dusty  particles  and  gas.  Because  of  lack  of  understanding  of  the 
drag  of  the  partic]e  cloud  and  heat  conduction,  current  analytical 
and  numerical  solutions  all  use  the  drag  and  heat  conduction  of  an 
isolated  rigid  sphere  at  steady  motion: 


JUf  ^ 

dt 

ili- 

dt 


4  D 

6KNu 

CD’p, 


—  (u  —  Uf)\u  —  Up\ 

Pf 


(T~Tp)  . 


(16) 

(17) 


The  parameters  in  the  relaxation  region  require  numerical  analysis. 
The  structure  of  shock  waves  in  dusty  gas,  according  to  the  experi¬ 
mental  conditions  of  experiment  No.  201,  is  shown  in  Fig.  1. 


When  the  coming-flow  velocity  UQ=aQ,  then  I,  and  the 

parameters  of  the  frozen  region  are  same  as  that  of  the  wave  front. 
Because  «« < ««.  and  —  the  parameters  of  equilibrium  are 

different  from  that  of  the  wave  front.  Under  this  situation,  disperse 
shock  waves  without  an  interruption  place  can  occur. 


Fig.  1.  Calculated  shock  wave 
structure  of  dusty  gas. 

(*/„  -  1.30,  7  -  0.17,  a  -  80;jm). 


t 


50  kpa ,  0.5  ms  per  div 

Fig.  2.  Typical  recorded  curves 

of  pressure.  Upper:  Without  dust, 

M  =1.44.  Lower  :  7— 0.17,  D  — sOfim,  M,,  -  1,24. 
s 


4.  Results  of  Observation 


As  described  before,  the  changes  of  shock  wave  velocity,  pressure 
and  dust  density  were  measured  in  the  experiments.  Because  the 
dusty  particles  were  gradually  absorbed  on  the  observation  window 
when  the  shock  waves  propagate  through  it,  the  measured  result  of 
particle  density  tended  to  be  a  little  higher.  For  the  particles 
under  free  fall  motion,  the  adsorption  was  not  serious  and  the 
initial  value  of  the  dust  density  is  quite  accurate.  We  now  give  a 
discussion  on  the  experimental  data  of  pressure,  combining  the 
shock  wave  velocity  and  the  initial  dust  density.  The  typical 
experimental  pressure  curves  are  shown  in  Fig.  2. 


1.  Pressure  in  the  frozen  and  equilibrium  regions.  According 
to  the  result  of  analysis,  the  pressure  in  the  frozen  region  Pi  is 
independent  of  the  loading  ratio,  size,  and  shape  of  the  dusty 
particles.  Fig.  3  is  the  experimental  result  with  »j  —  0—0.4,  A/,, — 
1.04—1.5,  D  —  SOfxm  and  55^m  .  The  curve  in  Fig.  3  was  calculated 
according  to  formula  (18).  From  this  figure  we  can  see  the  experi¬ 
mental  result  verifies  the  theoretical  prediction 


(Pi-P,) 

Po 


2r 

r  ■+■  1 


(  M!,  -  1  ]  . 


(18) 


Fig.  3.  The  relation  between  the  pressure  of  frozen  region  and  Mach 
number  of  shock  wave. 

Key:  (1).  without  dust. 
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Fig.  4a.  Calculated  results 
with  different  resistance 
coefficients . 


Fig.  4b.  Ratios  of  calculated 
values  to  experimentally  measured 
values . 

Key:  (1).  Experimental  curve. 


The  reason  for  the  disagreement  of  calculation  with  Stokes 
drag  from  the  experimental  result  is  obvious:  because  the  Reynolds 
number  under  the  above  condition  is  much  larger  than  one.  Although 
the  standard  curve  of  the  resistance  formula  can  represent  the  drag 
of  isolated  sphere  steady  motion,  it  can  not  reflect  the  influence 
of  the  disturbance  of  multi-particle  tail  flow,  particle  shape, 
surface  stature,  acceleration,  etc.  Because  the  measured  result  is 
far  from  the  calculated  result,  the  characteristics  of  the  relatation 
region  need  further  study. 

3.  Structure  of  disperse  shock  waves.  According  to  the  above 
analysis,  when  the  propagation  velocity  of  shock  waves  equals  the 
sound  velocity  of  the  gas,  disperse  shock  waves  without  interrupted 
wave  front  plane  will  appear.  Because  the  intensity  of  the  shock 
wave  is  weak  and  there  is  no  interrupted  plane,  it  is  quaite  difficult 
to  generate  and  observe  such  shock  waves.  We  recorded  the  pressure 
curve  of  such  shock  waves  in  the  experiment  of  No.  219  (Fig.  5). 

The  condition  of  this  experiment  was  as  follows:  driving 
gas  pressure  equals  0.54  MPa,  driven  gas  pressure  equals  0.1  MPa, 
dust  particle  diameter  equals  20-30  initial  dust  concentration 

equals  0.3  mg/cm^,  room  temperature  equals  27°C.  The  propagation 
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velocity  of  shock  waves  was  obtained  by  comparing  the  output  of  two 
pressure  sensors  which  are  located  a  certain  distance  apart.  The 
propagation  velocity  measured  was  355  m/s. 

4.  The  damping  effect  on  the  disturbance  of  dust  particles. 

The  upper  curve  in  Fig.  2  is  the  pressure  curve  without  dust  and  the 
lower  curve  is  the  pressure  with  dust,  and  all  other  experimental 
conditions  were  the  same.  We  can  see  that  the  fluctuation  of  the 
upper  curve  is  obviously  larger  than  that  of  the  lower  curve. 

Usually  such  fluctuation  is  caused  by  turbulent  flow  fluctuation  and 
the  unevenness  of  the  tube  wall.  The  dust  particles  in  gas  have 

an  obvious  effect  to  depress  such  fluctuation. 

5.  Driving  characteristics  of  dusty  gas  shock  wave  tube. 

Recently  different  studies  and  applications  using  dusty  gas  shock 
wave  tubes  are  increasing  gradually.  In  the  design  of  such  shock  wave 
tubes,  obtaining  the  required  driving  force  is  one  of  the  main 
concerned  problems. 

Only  if  the  propagation  distance  of  the  shock  wave,  that  is,  the 
distance  between  the  shock  waves  and  the  boundary  plane,  is  larger 
than  relaxation  length,  then  the  pre: sure  and  velocity  in  the 
boundary  plane  will  equal  the  pressure  and  velocity  in  the  equilibrium 
region.  The  calculation  of  structural  parameters  of  a  dusty  gas 
shock  wave  tube  can  be  carried  out  by  using  the  calculation  formulas 
of  the  classic  shock  wave  tube,  by  replacing  the  pure  gas  parameters 
with  the  thermal  dynamic  parameters  of  dusty  gas. 

Figure  6  gives  the  calculated  results  and  experimental  results 
of  the  variable  cross-section  shock  wave  tube  used  in  our  experiments. 
The  data  in  this  figure  are  the  shock  wave  velocities  in  the 
observation  region.  From  this  figure  we  know  the  calculated 
results  using  the  above  model  agree  with  the  experimental  results. 


25  kpa,  0.5  ms/division 

Fig.  5.  Pressure  curve  of 
disperse  shock  wave. 


Fig.  6.  Relation  of  particle 
loading  ratio  and  shock  wave  Mach 
number . 

Key:  (1).  Calculated  values 

without  dusty  particles.  (2). 
Measured  values  without  dusty 
particles.  (3) .  Measured  values 
of  dusty  gas. 


5.  Conclusion 


I 


1.  The  parameters  of  the  frozen  and  equilibrium  regions  can  be 
correctly  predicted  by  the  theoretical  analysis. 

2.  For  the  relaxation  region,  we  choose  the  drag  and  heat 
conductivity  of  an  isolated  rigid  sphere  at  steady  motion.  The 
tendency  of  change  is  same  as  that  of  the  experimental  result, 
but  the  quantitative  difference  if  very  large.  Therefore  further 
studies  are  needed. 

3.  When  shock  waves  propagate  with  the  sound  velocity  of  gas 
in  a  wave  front,  disperse  shock  waves  without  interrupted  plane 
were  observed.  This  verifies  the  theoretical  prediction. 

4.  The  dusty  particles  have  the  disturbance  effect  in  the 
damping  gas  flow. 

5.  At  the  location  far  away  from  the  diaphragm  the  shock  wave 
intensity  can  be  predicted  by  using  the  classic  shock  wave  formulas. 
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with  the  effective  gas  thermal  dynamic  parameters. 


The  authors  wish  to  thank  Mr.  Li  Zhongfa  and  Mr.  Pu  Yi  Kang  for 
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